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Drag reduction by polymers is bounded between two universal asymptotes, the von-Karman log- 
law of the law and the Maximum Drag Reduction (MDR) asymptote. It is theoretically understood 
why the MDR asymptote is universal, independent of whether the polymers are flexible or rodlike. 
The cross-over behavior from the Newtonian von-Karman log-law to the MDR is however not uni- 
versal, showing different characteristics for flexible and rodlike polymers. In this paper we provide 
a theory for this cross-over phenomenology. 



I. INTRODUCTION 

The phenomenon known as "draff reduction" by poly- 
mers in turbulent channel flows [ij, |2j is conveniently dis- 
cussed in channel geometry for fixed pressure heads, such 
that the reduction in the drag is manifested as an in- 
creased mean velocity. For the sake of comparison be- 
tween different fluids it is convenient to choose normal- 
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FIG. 1: Typical velocity profiles taken from Q. In dashed 
line we noted the von-Karman law while the MDR © 
is shown as the continuous black line. In all cases the mean 
velocity follows the same viscous behavior for y + < 10. After 
that the scenario is different for flexible and rodlike poly- 
mers. The typical behavior for the former is presented by the 
open triangles, which follow the MDR up to a crossover point 
that depends on the concentration of the polymer and on the 
value of Re. The rodlike behavior is exemplified by the solid 
triangles and the open squares; the mean velocity profiles ap- 
pear to interpolate smoothly between the two asymptotes as 
a function of the concentration of the rodlike polymer. 



ized coordinates. Denote the mean pressure (per unit 
density) gradient p' = —dp/dx where x, y and z are the 
lengthwise, wall-normal and spanwise directions respec- 
tively. The length and width of the channel are usually 
taken much larger than the mid-channel height L, mak- 
ing the latter a natural re-scaling length for the introduc- 
tion of dimensionless (similarity) variables, also known as 
"wall units" Q. Thus the Reynolds number Re, the nor- 
malized distance from the wall y + and the normalized 
mean velocity V + (y + ) (which is in the x direction with 
a dependence on y only) are defined by 

Re = L^/Pl/vq , y + = yRe/L , V + = V/^/p^L , (1) 

where vq is the kinematic viscosity of the neat fluid. For 
Newtonian fluids the profile of the mean velocity V + {y + ) 
is universal, starting with the viscous sub-layer in which 
V + (y + ) = y + and then, at y + somewhere between 6 and 
12, the profile crosses over to the universal von-Karman 
log-law of the wall (cf. Fig. ^) 

V + {y + ) = n- l \ny+ + B . (2) 

Upon the addition of small concentrations of polymers, 
the drag is reduced and for the same value of p' one 
finds an increase in V + (y + ). This phenomenon exhibits 
both universal and non-universal aspects. The universal 
aspect is the Maximum Drag Reduction (MDR) asymp- 
tote, which is the largest attainable profile V + (y + ). This 
was determined experimentally by Virk who found |2| 

V + {y + ) = — ln(e/c v y+) for y+ > 12 . (3) 

While n^ 1 2.5, k^ 1 w 12, leading to a significantly 
larger mean velocity at the MDR as compared with von- 
Karman's log-law. The MDR appears independent of the 
nature of the polymer (for example whether it is flexible 
or rodlike), of the length of the polymer and of the con- 
centration. On the other hand, the way that the system 
attains the MDR is not universal, and it depends on all 
of the above. In particular, it appears that flexible and 
rodlike polymers attain the MDR, as a function of the 
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concentration, in qualitatively different ways (see for ex- 
ample |2|). The experimental information is quite scant, 
but available data indicate different scenarios for the two 
types of polymers. The data in Fig. 2] indicate for large 
values of Re the mean velocity profile with flexible poly- 
mers [polyacrylamide (PAA)] follows the MDR until a 
point of cross over back to the "Newtonian plug" , where 
it becomes parallel to von-Karman's log-law. Increasing 
the concentration results in following the MDR further 
until a higher cross over point is attained back to the 
Newtonian plug 2] . On the other hand, for rodlike poly- 
mers [sodium carboxymethylcellulose (CMC) and sodium 
carboxymethylccllulosc/xanthan gum blend (CMC/XG)] 
the data shown in Fig. ^ indicate a different scenario. 
Contrary to flexible polymers, here, as a function of the 
concentration, one finds mean velocity profiles that inter- 
polate between the two asymptotes J2J and , reaching 
the MDR only for large concentrations. A difference in 
the behavior of flexible and rodlike polymers when plot- 
ting the drag as a function of Reynolds number was also 
reported by Virk and coworkers, see JjJ. 

The universal MDR was fully explained in recent work, 
and its parameters derived, by showing that it is an edge 
turbulent solution in a channel @. In other words, if 
one tried to reduce the drag further (or further increase 
the profile V + (y + )), one would lose the turbulent flow in 
favor of a laminar solution. This is the reason for the uni- 
versality of the MDR and its insensitivity to the nature of 
the polymer. The aim of this paper is to derive quantita- 
tively the non-universal scenarios of attaining the MDR 
by flexible and rodlike polymers. We will limit our at- 
tention to the case of high Re, and consider the profiles 
V + {y + ) for different concentrations, with the aim of ex- 
plaining the phenomenology displayed in Fig. ^ The pa- 
per ends with a definite prediction of different crossover 
behavior in the mean velocity profiles for flexible and rod- 
like polymers, in agreement with the indication of Fig. ^ 

In Sect. 2 we summarize the available theory of drag 
reduction by flexible and rodlike polymers, based on the 
balance equations for energy and momentum in the tur- 
bulent boundary layer. We derive in this section the 
equations satisfied by the y dependence of the mean shear 
(in wall units), from which the velocity profiles are ob- 
tained by integration. Sect. 3 presents the results to- 
gether with a summary and a discussion. 



II. THEORY OF DRAG REDUCTION 



A. The polymeric stress tensor 



The extra stress tensor er is due to the interaction be- 
tween the polymers and the fluid. Within the FENE-P 
model for flexible polymers we have Q 



Cab ~ V p J p lZ a b , 



(5) 



where v p is the polymeric contribution to the viscosity 
at zero shear, and 7 P is the inverse relaxation time of 
the stretched polymer. In writing this expression one 
adopts the standard simplification of a single relaxation 
time. The conformation tensor 1Z is obtained from the 
normalized end-to-end distance vector of the polymer 
p = p/pmax, averaged over the conformation of the poly- 
mers, 



K= pp 



(6) 



For rodlike polymer p becomes a unit vector, and there 
is no coil-stretching transition. Accordingly the stress 
tensor assumes a different form m : 



Cab — ftVplZaglZijSi 



(7) 



where SVj is the strain experienced by the polymer = 
dUi/drj. As explained above, the difference in form of 
the stress tensor is immaterial for the universal form of 
the MDR 0, El , a phenomenon that was called "additive 
equivalence" by Virk. We will see that this difference 
translates however to a very different scenario for the 
attainment of the MDR. 



B. The balance equations 

The phenomenon of drag reduction can be understood 
on the basis of the balance equations of the mechanical 
momentum and turbulent energy |lCj . These are derived 
on the basis of the Reynolds decomposition 



Ui(r,t) = V(y)5 ix +Ui(r,t) 



(8) 



S l3 (r , t) = S(y)S lx S jy + s tJ (r , i) , S(y) = (.9) 

Writing these equations we use the fact that in a turbu- 
lent channel flowp' is constant, and due to the symmetry 
all the other mean quantities depend on y only. In ad- 
dition to the mean shear S(y), we need to to introduce 



the mean turbulent kinetic energy K 



2 ) /2 and the 



Reynolds stress W = —(u x u y ). The momentum balance 
equation is obtained by averaging Eq. (0J and integrat- 
ing with respect to y, ending up with the exact equation: 



In the presence of a small concentration of polymers 
the Navier-Stokes equations for the fluid velocity U(r, t) 
gain an additional stress tensor: 



dU 

dt 



U-VU 

V • u 



-Vp + u V 2 U + V • <t , (4) 







i> a S + W + (a xy ) = p'(L-y) . 



(10) 



In wall-units Eq. (|10|l can be written in a more elegant 
form: 

S+ + W+ + (<r+ ) = (1 - y + /Re) (11) 
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where S+ = v S/(p'L), W+ = W/(p'L), and a+ = 
Oij/{p'L). When Re is very large and for y + not too 
large we neglect the second term on the RHS, approxi- 
mating the RHS as unity. 

The balance equation for the turbulent kinetic energy 
is calculated by taking the dot product of the fluctuation 
part of Eq. @ with u: 



d 

WS = —(UyU 2 + UyP ~ (TiyUi) 



V()(Si 



(12) 

Also this equation is exact. We simplify it by noting 
that the first term on the RHS involves the spatial flux 
of turbulent energy which is known to be negligible in 
the log-layer. The second term represents the dissipation 
which is modelled (in wall- units) as [Toj : 



( s ij s ij) 



a (K+) 3 / 2 

/;+ 11+ 



(13) 



where = h / os i j/(p'L) and K + = K/(p'L). Therefore, 
the energy balance equation takes on the form 



W + S A 



V + 



(14) 



Finally, we quote the experimental fact that in the log- 
layer W + and K + are proportional to each other: 



z 2 = w + 



(15) 



Experimentally, it was found that c ~ 0.5 in the Newto- 
nian case, and c w 0.25 in the MDR. 



C. Effect of the polymers 

In the Newtonian case, when a = and c = 0.5, the 
three equations (fill) . I|14|) and (|15fl are sufficient for de- 
termining the three unknowns S + , K + and W + . The 
best fit to Newtonian experiments and simulations are 
obtained using the values a = 3 and b = 0.321, and we 
are going to use these values throughout, also in the vis- 
coelastic cases discussed below. In the presence of poly- 
mers, however, we have to consider the terms introduced 
by the polymers to these equations. The necessary the- 
ory was presented in [flLllflj. with the final results relating 
these terms to the yy component of the mean conforma- 
tion tensor R = (TV : . 



(<Jxy) 
((JijSij) 



K 



(16) 



Using these theoretical results the momentum and energy 
equations, in the limit of large Re, are reduced to 



S+ + W+ + vR vv S + = 1 



(17) 



and 



K 



y + y + 



(18) 



Here v = v p /vq. To proceed we need to relate R yy to 
the other variables. The necessary theory is presented in 
E3 El EH leading to an important difference between 
the flexible and the rodlike polymers: 



R 



ny 



v K 

(flexible) , 

Sy 



K 



(rodlike) 



yy S 2 y 2 
Substituting into Eqs. l(T7Jl and (fT5f) we have 
u cS S+ +W+ = 1, 

and 



W + S+ = K + v eS (^) 2 +b 



(K+) 



3/2 



with the "effective viscosity" 




(flexible) , 
(rodlike) . 



(19) 
(20) 

(21) 

(22) 

(23) 
(24) 



In the next section we will show that although the ef- 
fective viscosities take on different forms for the flexible 
and rodlike cases, in fact they both depend linearly on y + 
whenever we have a log layer with S + oc 1 /y + . The rea- 
son is that K + turns out to be proportional to (y + ) 2 and 
to y + for the flexible and the rodlike cases respectively. 
Accordingly, we will write in both cases 



v eS = 1 + a[y + - A(a)] 



(25) 



with A(a) being the width of the viscous sub-layer, and 
its dependence on the slope of the effective viscosity a 
needs to be determined. It is natural to present A(a) in 
terms of a dimensionless scaling function f(x), 



A(a) = 5 + f(aS+) , 



(26) 



where S + w 6 is the width of the Newtonian viscous 
boundary layer. In the Newtonian limit a — > 0, p e ff ~ y 1 
and A — > <5 + , hence we have /(0) = 1. In |6j it was 
shown that the balance equations l|21|) and (1221 (with the 
prescribed form of the effective viscosity profile) have an 
non-trivial symmetry that leaves them invariant under 
rescaling of the wall units. This symmetry dictates the 
function A(a) in the form 



A(a) 



1 - a5+ 



(27) 
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D. Closing the equations 
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To complete the model, we have to specify the value 
of c in Eq. 115|) . This parameter becomes naturally a 
function of a. We can find its a-dependence by identi- 
fying the width of the viscous sub-layer A with a/c(a). 
This stems from the fact that the balance equations can- 
not support a turbulent solution for y + < a/c(a). This 
means that 



a/c(a) = A(a) . 



(28) 



Combining then Eqs. (HSJ), (H3 and and putting 

5 + = 6, wc can solve and find 



C(>eff) 



a a(2 - v eS ) 
6 



2yH 



a 

12 



/ 1 12^ cff | 36(2-^ cff )2 i 



IT 



(y 



+ \2 



(29) 



To summarize note that Eqs. (|23|l and (|24|l can be 
written as 



A\S + y + ) 



+ \2 



(30) 
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FIG. 2: The mean velocity profiles for flexible polymer ad- 
ditives with v = 1,5,10,20,50,100 and 500 from below to 
above. Note that the profile follows the MDR until it crosses 
over back to the Newtonian plug. 
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with 



A 2 = 
A 2 = 



v e s - 1 

V 

V 



(flexible) 
(rodlike) 



(31) 
(32) 



Using Eqs. (|15fl and (|30|l . we can rewrite Eqs. (|21|l and 
(1221) as two equations for the two variables z^ e ff and 5 I+ : 



v eS S + + c 2 A 2 (S + y + f = l 



and 



2 Q+ 



v eS (—y + bAs^ 



Equation (|31Jl implies 



S+ 



(y+) 2 {c 2 -bA) 



(33) 



(34) 



(35) 



Substituting Eq. ^ into Eq. ^ gives an equation for 
vl S (^?(c 2 - bA) + c 2 A 2 v 2 a (^) 2 = (c 2 - bAf (36) 

y + v + 

Finally, we can solve Eq. (|3l)|) to get v e s{y + ) for dif- 
ferent values of v. Then we can obtain S + and K + using 
Eqs. (|33|l and lj5U|) respectively. Integrating S + over y + , 
we get also V + (y + ). 
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FIG. 3: The mean velocity profile for rodlike polymer addi- 
tives with v = 1,5,10,20,50,100,500,1000,5000 and 10000 
from below to above. Note the typical behavior expected for 
rodlike polyemrs, i.e. that the profile diverges from the von 
Karman log-law, reaching the MDR only asymptotically. 



III. RESULTS 

The results of the numerical solutions of the equations 
are shown in Figs. J2HBland0] 

In Figs. |21 and 13 we present the mean velocity profile as 
a function of the distance from the wall, for flexible and 
rodlikc polymers respectively. The main result of this 
paper is seen in the difference between these profiles as a 
function of the polymer concentration. While the flexible 
polymer case exhibits the feature 0, that the velocity 
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FIG. 4: The turbulent kinetic energy profile K + {y + ) for the 
flexible (solid line) and rodlike (dashed line) cases respec- 
tively. We note the quadratic and linear dependence respec- 
tively, as anticipated in the text. 



values of v. This result is in agreement with the exper- 
imental finding in 0, ^] that the flexible polymer is a 
better drag reducer than the rodlike analogue. 

We should note that the higher efficacy of flexible poly- 
mers cannot be easily related to their elongational vis- 
cosity as measured in laminar flows. In some studies 
(l2l fill. IT^ I it was proposed that there is a correlation 
between the elongational viscosity measured in laminar 
flows and the drag reduction measured in turbulent flows. 
We find here that flexible polymers do better in turbu- 
lent flows due to their contribution to the effective shear 
viscosity, and their improved capability in drag reduction 
stems simply from their ability to stretch, something that 
rodlike polymers cannot do. 

Finally we recall that the derivation of our equations 
relied on the fact that the turbulent kinetic energy K + 
depends linearly and quadratically on the distance from 
the wall for flexible and rodlike polymers respectively. It 
is important to check that the resulting equations confirm 
this expectation self-consistently. Indeed, in Fig. 0|we 
present the solution for this quantity in the two cases, 
and find that the expectation is fully realized. 



profile adheres to the MDR until a crossover to the New- 
tonian plug is realized, the rodlike case presents a "di- 
verging" of profiles which only asymptotically reach the 
MDR. We also notice that the flexible polymer matches 
the MDR faithfully for relatively low values of v, whereas 
the rodlike case attains the MDR only for much higher 
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